function [dPsiEta] = dPsiPNbydEta(f, Params)
% 
% DPSIPNBYETA - compute the derivative of the 3.5PN phase with respect
% to the symmetric mass ratio eta. 
% 
% usage: [dPsiEta] = dPsiPNbydEta(f, Params)
%
% P. Ajith, 05.04.2008
%
% $Id: dPsiPNbydEta.m 59 2010-01-20 21:41:01Z anand.sengupta $

    MSOLAR_TIME = 4.92579497077314e-06;

    m = Params.totalMass.*MSOLAR_TIME;
    eta = Params.eta;

dPsiEta = 3085./(3072.*power(f.*m,0.3333333333333333).*power(pi,0.3333333333333333)) + ...
  (76055.*power(f.*m,0.3333333333333333).*power(pi,0.3333333333333333))./...
   73728. - (74045.*power(f.*m,0.6666666666666666).*...
     power(pi,1.6666666666666667))./32256. + ...
  (9.551407715020614.*power(f.*m,0.3333333333333333))./power(eta,2) - ...
  3./(128.*power(f.*m,1.6666666666666667).*power(pi,1.6666666666666667).*...
     power(eta,2)) - 3715./(32256.*f.*m.*pi.*power(eta,2)) - ...
  15293365./(2.1676032e7.*power(f.*m,0.3333333333333333).*...
     power(pi,0.3333333333333333).*power(eta,2)) + ...
  (3.*power(pi,0.3333333333333333))./...
   (8.*power(f.*m,0.6666666666666666).*power(eta,2)) - ...
  (38645.*pi)./(32256.*power(eta,2)) - ...
  (77096675.*power(f.*m,0.6666666666666666).*power(pi,1.6666666666666667))./...
   (1.0838016e7.*power(eta,2)) - ...
  (127825.*power(f.*m,0.3333333333333333).*power(pi,0.3333333333333333).*...
     eta)./27648. - (38645.*pi.*log(6))./(21504.*power(eta,2)) + ...
  (107.*power(f.*m,0.3333333333333333).*power(pi,0.3333333333333333).*...
     log(power(f.*m,0.3333333333333333).*power(pi,0.3333333333333333)))./...
   (14.*power(eta,2)) - (38645.*pi.*...
     log(power(f.*m,0.3333333333333333).*power(pi,0.3333333333333333)))./...
   (10752.*power(eta,2));
